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Abstract — This paper presents a novel two-way decode-and- 
forward (DF) relay strategy for Orthogonal Frequency Division 
Multiplexing (OFDM) relay networks. This DF relay strategy em- 
ploys multi-subcarrier joint channel coding to leverage frequency 
selective fading, and thus can achieve a higher data rate than 
the conventional per-subcarrier DF relay strategies. We further 
propose a low-complexity, optimal power allocation strategy to 
maximize the data rate of the proposed relay strategy. Simulation 
results suggest that our strategy obtains a substantial gain over 
the per-subcarrier DF relay strategies, and also outperforms 
the amplify-and-forward (AF) relay strategy in a wide signal- 
to-noise-ratio (SNR) region. 

I. Introduction 

In recent years, relaying has emerged as a powerful tech- 
nique to improve the coverage and throughput of wireless 
networks. Compared with the traditional one-way relaying, 
two-way relaying provides better spectral efficiency, where 
two terminal nodes employs an intermediate relay node to 
exchange information simultaneously 0,0. 

Orthogonal Frequency Division Multiplexing (OFDM) is 
an essential broadband transmission technique to improve the 
spectral efficiency of wireless networks. A combination of 
OFDM and relaying techniques has been advocated by many 
industry standardization groups of next generation wireless 
networks, such as IEEE 802.16m and 3GPP's LET- Advanced. 

In one-way OFDM relay networks, multi-subcarrier joint 
decode-and-forward (DF) relaying was studied in 0-0, 
which can achieve higher data rate than per-subcarrier DF 
relaying. For two-way OFDM relay networks, the amplify- 
and-forward (AF) relay strategies were commonly adopted 
[6|-|8|. However, their performance is quite poor in the low 
signal-to-noise-ratio (SNR) region due to the amplified noises. 
The per-subcarrier DF relay strategies were considered in ||9l- 
[11 1, which are essentially simple accumulations of narrow- 
band two-way DF relaying over the individual subcarriers. 
Unfortunately, these strategies suffer from rate losses due to 
channel mismatching. 
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In this paper, we propose a novel multi-subcarrier DF 
relay strategy for two-way OFDM relay networks. By per- 
forming channel coding across subcarriers, this strategy can 
exploit frequency selective fading, and achieve higher data 
rate than the per-subcarrier DF relay strategy in J9j- We 
further formulate a power allocation problem to maximize 
the exchange rate, which is defined as the maximal data 
rate can be simultaneously achieved in both directions. An 
efficient dual decomposition algorithm is proposed to resolve 
this problem, which has a linear complexity with respect to 
the number of subcarriers. Simulation results show that the 
proposed multi-subcarrier DF relay strategy outperforms not 
only the conventional per-subcarrier DF relay strategy, but also 
the AF relay strategy in a wide SNR region. 

II. System Description 

We consider the two-way OFDM relay network shown in 
Fig. [T] two terminal nodes T\ and exchange information 
via an intermediate relay node Tr. Assume that each node 
has a single antenna and operates in a half-duplex mode, i.e., 
transmitting and receiving in orthogonal time slots 0, 0- 
All the nodes employ OFDM air interface with the same N 
subcarriers. 
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(a) The MA phase. 
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(b) The BC phase. 
Fig. 1. System model of two-way OFDM relay network. 



The DF relay procedure comprises of a multiple-access 
(MA) phase and a broadcast (BC) phase without direct trans- 
missions, as shown in Fig. Q] We set X{ = (Xn, . . . , Xyy) 
and Yi = (Y a ,. ■ ■ , Y iN ) for % = 1, 2, i?, where Xj„ and F in 
denote the normalized transmitted and received signal in the 
nth subcarrier at T,;. In the MA phase, T\ and T 2 transmit 
X\ and X2 simultaneously to the relay node, and the relay 
Tr performs multi-user detection to fully decode X\ and X2 
from the received Y R ; in the BC phase, the relay T R broadcasts 
X R = f(X 1 ,X 2 )toT ia ndT 2 . 

Specifically, in the nth subcarrier, h ln and h 2n denote the 
channel coefficients from T\ and T 2 to T R , respectively, h\ n 
and /i2n denote the channel coefficients from T R to T\ and 
T2, respectively. Thus, the received signals Yi n 's in the nth 
subcarrier at T,'s are given by 

^Rn = \J Plnh\ n Xi n + \J P 2n h 2n X 2n + Z Rn , (1) 
Yin = yPRnhlnXjin + Z\ n , (2) 

Y 2n = \J PR n h 2n Xji n + Z 2n , (3) 

where Pj, n denotes the transmit power in the nth subcarrier 
at Ti, and Z^ n denotes independent complex additive white 
Gaussian noises with zero mean and unit variance, i.e., Z$ n ~ 
CAf (0,1), for i = 1,2, R. Therefore, P in essentially denotes 
the corresponding transmit SNR. 

We assume that p, 6 (0, 1) denotes the fixed proportion of 
time slot allocated to the MA phase, and all the terminals 
are subject to separate power constraints X^Li-Pin < Pimax 
(i = 1,2, R), where Pi max denotes the maximum available 
power for Ti. 

III. A NOVEL DF RELAY STRATEGY FOR TWO-WAY 
OFDM RELAY NETWORKS 

For two-way OFDM relay networks, the conventional DF 
relay strategies simply applied the narrow-band DF technique 
over each subcarrier independently, and the overall throughput 
was the sum rate of all the subcarriers (9), ifTTI . However, these 
strategies suffer from rate losses due to channel mismatching. 

In this section, we propose a novel multi- sub carrier DF 
relay strategy, which performs channel coding across sub- 
carriers to leverage frequency selective fading, and achieve 
higher data rate. An achievable rate region for this strategy is 
provided in the following theorem. 

Theorem 1. Any rate pair (R\ 2 , R 2 \) satisfying the follow- 
ing inequalities is achievable for the two-way OFDM relay 
network given by ([T)-(|3j: 

N 

R12 <A*5^1og a (l + |^in| 2 i\n), (4a) 

n=l 

N 

Rl2 < (1 - A») 1 °8'2( 1 + \~ h 2n\ 2 PRn), (4b) 

n=l 

N 

R21 <V^\og 2 (l + \h 2n \ 2 P 2n ), (4c) 

71=1 



N 

R21 < (l-u)^log 2 (l + |/H„| 2 PR n ), (4d) 

n=l 

N 

R12 + R21 < log 2 (l + \h ln \ 2 P ln + \h 2n \ 2 P 2n ), (4e) 

n=l 

N 

Pin < Pimax, i = l,2,R, (4f) 

n=l 

Pin>0, VnGAT, i=l,2,R. (4g) 

where Af = {1, . . . , N}, R\ 2 and R 2 \ denote the achievable 
data rates from T\ to T 2 and from T 2 to T\, respectively. 

Proof: An achievable rate region of a DF relay strategy 
for the discrete-memoryless two-way relay network has been 
given by the set of (R\ 2 , R21) rate pairs satisfying lfl2l 

R12 < min{/i/(X i; Y R \X 2 ), (1 - p,)I(X R ; Y 2 )} , (5a) 
R21 < min{nI(X 2 ; Y R \X X ), (1 - n)I(X R ; Y x )} , (5b) 
Ri2 + R2i<fiI(X 1 ,X 2 ;Y R ). (5c) 

Each mutual information item in (f5ab-(f5cl> corresponds to the 
achievable rate of a parallel point-to-point channel. Similar 
to the idea in the proof of [3 Theorem 1], we choose the 
input signals for each subcarrier to be independent Gaussian 
distributed with unit variance, i.e., Xi n ^CN(Q, 1). Thus, each 
mutual information item in (fSat-dScli is replaced by the sum 
of N logarithmic rate items decided by (HJ-© with separate 
power constraints. The proof is complete. ■ 
Remark 1: The key idea of this multi-subcarrier two-way 
DF relay strategy is introducing channel coding across sub- 
carriers to fully exploit frequency selective fading. Implicitly, 
the information transmitted over one subcarrier in the MA 
phase may be forwarded over some other subcarriers in the 
BC phase. By this, the problem from mismatching of wireless 
channels over subcarriers is resolved. The achievable rate 
region of Theorem [T] is no smaller than that achieved by the 
per-subcarrier two-way DF relaying, which is the set of rate 
pairs satisfying J5] 

N 

R\2 < y^min{/ilog 2 (l + \h ln \ 2 P ln ) , 
n= \l-[i)log 2 (l + \h 2n \ 2 P Rn )}, 

N 

R21 < ^min{ulog 2 (l + \h 2n \ 2 P 2n ) , 
n= \l-[i)log 2 (l + \h ln \ 2 P Rn )}, 

N 

R12 + R21 <//^log 2 (l + |/ii n\ 2 Pln + \h 2n \ 2 P 2n ). 

n=l 

Therefore, multi-subcarrier two-way relay channel is not 
a simple linear combination of multiple narrow-band single- 
subcarrier two-way relay subchannels. Similar observations 
have been found for one-way parallel relay networks Q. 

IV. Optimal Power Allocation 

In this section, we investigate the largest achievable sym- 
metric exchange data rate of our proposed DF relay stratey, 
which can be approached by an optimal power allocation. 



A. Problem Formulation 

By optimizing the power allocation strategy (Pi, P 2 , Pr), 
our objective is to maximize the exchange rate Rx = 
min{i?i2, -R21}, which is defined as the data rate can be 
achieved simultaneously in both directions, where P; = 
[Pji, Pi2, . . • , Pjjv] T denotes the power allocation vector at Ti, 
for i = l,2,R. This can be expressed as the following convex 
optimization problem: 



max i?x 

P u P 2 ,Pr,Rx 



N 



(7a) 
(7b) 



s.t. R x < //^log 2 (l + \h ln \ 2 P ln ), 
n=l 

N 

Rx < (l-^E^^ + l^npPftn), (7c) 

71=1 

JV 

i?x < fi^2log 2 (l + \h 2n \ 2 P 2n ), (7d) 

n=l 

AT 

i?X < (1 - /i) J] lo S2 (l + l^ln l^fln) , (7e) 



N 



Rx<^J2 lo g2(! + \hm\ 2 Pm + \h 2n \ 2 P 2n ), (7f) 



JV 



^ ^ Pn Pmax; ^ — 1? 2, .R, 
n=l 

Pin > 0, V n e AT, i = l,2,i?. 



(7g) 
(7h) 



It is readily observed that in the problem (O, Pi and P2 are 
only related to the constraints (l7bl (fTdl i ( fTfT i. while Pr is only 
related to the constraints (TTcb ( TTeb . This observation helps to 
decompose our original power allocation problem © into the 
following two subproblems: 



max Rma 
p u p 2 ,Rua 



(8) 



JV 



s.t. Rma < ^E 10 ^ 1 + l^m| 2 -Pin), 

71=1 

AT 

Rma < \i lo §2 (1 + |ftan| 2 ftn) » 

71 = 1 

iV 

Rma < ^ ^log 2 (l + |/ii„| 2 Pi„ + |^ 2 „| 2 P„) 



JV 



iV 



fin < flmax, Pn < P 



2maxi 



Pm>0,P 2 „>0, V n e Af. 



max i?sc 

Pr,Hbc 



(9) 



S.t. P BC < (1 -/i)E 1 °S2( 1 + l^ln| 2 PRn 
n=l 

JV 

Pbc < (1 -//) E ^(l + IM 2 Pr„ 

n=l 

JV 

^ ^Rn < fflraax, Pfln > 0, V n G A/". 



We can denote -R^A an£ l ^bc as me optimal values for 
the MA subproblem ([8]) and the BC subproblem (O, respec- 
tively. Eventually, the maximal practical exchange rate for our 
proposed DF strategy is given by R x = mm{R MA , P BC }. 



B. Proposed Dual Decomposition Algorithm 

The interior-point methods can be used to solve both of 
the convex optimization problems ([8]) and (O, however, they 
quickly become computationally intractable as N increases, 
because they have a complexity of 0(N 3 ) at least when 
solving the search direction in each iteration lT3ll . Therefore, 
we present a low-complexity dual decomposition algorithm for 
the subproblems ^ and (O, to efficiently obtain the optimal 
solution to (|7]l. Next, we will take the subproblem (0 as an 
example to illustrate this algorithm. 

Note that problem ([8]) is strictly feasible. Then, according to 
the Slater's condition lfT3l . it is equivalent with the following 
dual optimization problem: 



max < min C (Pi , P2 , Rma ,\ot) 



(10) 



where 



£(Pi,P 2 ,P M A,A,a) 

JV 

Rma - /i^log 2 (l + \h ln \ 2 P ln ) 



+ Ai 
+ A 2 



n=l 



Rma- LiJ2log 2 (l + \h 2n \ 2 P 2 n) 



+ A3 


-Rma - 




/ N 


+ OL\ 




JV r 


\7l=l 
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71=1 L 


aiPin - 
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71=1 

JV 



71=1 
Pin 



N 



■ OL 2 \ 



E^ 



p 2 , 



eeiPin + a 2 P 2n - /xAi log 2 (l + \h\ n \ Pi„) 



^A 2 log 2 (l + \h 2n \ 2 P 2n ) 

^log 2 (l + \h ln \ 2 P ln + \h 2n \ 2 P 2n ) 



+ (Ai + A 2 + A 3 - 1) Rma - aiPin 



a 2 P 2 



(11) 



is the partial Lagrangian of ([8]), and A = [Ai, A 2 , A3] 11 , ot — 
[ai,a 2 ] T are nonnegative dual variables associated with the 
three rate constraints and two power constraints, respectively. 

According to ( fTTT l. the inner minimization problem of ( TTOb 
can be decomposed as N independent per-subcarrier power 
allocation problems. Hence, the computational complexity 
for solving the inner problem is only linear with respect 
to N. In addition, the optimal (Pi n ,P 2n ) must satisfy the 
following Karush-Kuhn-Tucker (KKT) conditions for given 



dual variables (A, a) |fl3l : 



dC 



0P Ul " l 21n2(l + \h ln \ 2 P ln + \h 2 „\ 2 P2n) 
»\i\h ln \ 2 f>0 if Pl„ = 



dC 



= a 2 - 



ln2(l + |/ ll „| 2 P 1 „) 1 = if P ln >0 

^A 3 |/t 2 »| 2 

21n2(l+|/ lln | 2 P ln + |/i 2n | 2 P 2n ) 
A*A 2 |fr 2n | 2 f >0 if P 2 „. = 



,(12) 



In 2(1 + |/i 2 „| 2 P 2rl ) 1 = if P 2 „ > 



(13) 



Thus, it must belong to one of the following four cases: 
Case 1: P Xn > 0,P 2n > 0. Then the formulas Q2) 
and ([T3l hold with equality. It is hard to solve 021 and 
(fT3l l directly since they are both quadratic equations of two 
variables P\ n and P 2n . However, we can utilize an auxiliary 



variable defined as x — \hi n \ 2 Pi r 



l P 2n to simplify 



them. More specifically, from (fT2l and Q3) , one can obtain 
that 



\h\ n \ 2 Pln — 



2/xAi|/li„| 



\h 2n \ 2 P 2 . 



2\n2- ai -fiX 3 \h ln \ 2 /(l + x) 

2ti\ 2 \h 2n \ 2 

21n2-a 2 - M A 3 |/ ) , 2 „| 2 /(l + ^) 



1, (14) 
1. (15) 



Taking the sum of the above two equations, we obtain a 
cubic equation of x, which has closed-form solutions given 
by Cardano's Formula [14]. After deriving the positive root x 
of this cubic equation, we can easily obtain the optimal P\ n 
and P 2n from 04) and 05) . By this procedure, the quadratic 
equations 02) and (fT~3T > are solved analytically by converting 
to an equivalent cubic equation. Finally, we need to check 
whether P\ n and P 2n satisfy the conditions P\ n > 0, P 2n > 0. 

Case 2: P Xn > 0, P 2 „ = 0. Then the solutions to ([T2j and 
(fT3l > can be derived as 



Plr, 



y(2\i + A3) 
2 In 2 • ai 



1 



P 2 „ = 0. 



(16) 
(17) 



This case happens only if P\ n > and the KKT condition 
03), 21n2-a 2 > 2^\ 2 \h 2n \ 2 + m^j^; , is satisfied. 

Case 3: Pi„ = 0, P 2n > 0. Then the KKT conditions can 
be reformulated as 



Pin = 0, 

M (2A 2 + A 3 ) 



Pin — 



1 



21n2-a 2 \h 2n \ 



(18) 
(19) 



This case happens only if P 2n > and the KKT condition 



-, is satisfied. 



(Ell, 21n2-«i > 2 f i\ 1 \h ln \ 2 + 1+{h2ril 2 P2ri 

Case 4: Pi„ = 0, P 2n — 0. This is the default case when 
the above three cases do not happen. 

Then, we optimize the dual variables (A, a) for the 
outer maximization problem of ( fTOb . We redefine v = 



[Ai, A 2 , A3, ai, a 2 ] T . Further, considering the KKT condition 
for the optimal data rate Pma, we have 



dC 



OR 



= Ai + A 2 + A 3 - 1 = 0. 



(20) 



MA 



In view of that the objective function is not differentiable 
with respect to (A, a), we consider to update v using the 
subgradient method fl5l . 0161 . Specifically, in the fcth iteration, 
the subgradient method updates u k by 



v k+1 = [u k + s k r,{u k )}. 



(21) 



where [v]-p represents the orthogonal projection of v to the 
dual feasible set {y \ 1 T A = l,X,a >z 0} based on a finite 
algorithm in ATI . s k is the step size of the fcth iteration, and 
r)(v k ) is the subgradient of the outer problem of ( fT0T > at v k , 
which can be chosen as 



n{v k )= 



~^n=^g 2 {l+\h ln \ 2 P* n ) 
-MElll0g2(l + l^2„| 2 P 2 * n ) 

-f Eli log a (l + \hm\ 2 PL + \h 2n \ 2 P 2 * n 

EN p* _ p 
n=l In ^lmax 

,N 



E 



P* 

= 1 



P 2 , 



,(22) 



where P* n and P% n are the optimal solution of the inner 
minimization problem in the fcth iteration. It has been shown 
that the subgradient updates in ( f2Tb can converge to the 
optimal dual point v* as fc —> 00, provided that the step size 
s is chosen according to a diminishing step size rule Ifl6| . 

Let Cma,i(Pi, P%){i = 1,2,3) denote the the right-hand 
sides of three rate constraints in (|8), respectively, and thus we 
obtain the optimal PJ^A = m i n {CMA,i(Pi*, ^2)1 * = 1> 2, 3}. 

The proposed dual decomposition algorithms for the MA 
subproblem ([8]) are summarized in Algorithm Q] Similarly, the 
BC subproblem (O can be solved with the same techniques. 
Their complexity grow in the order of O(N), which are much 
lower than the classic convex optimization software package 
based on interior-point methods. Therefore, our proposed 
algorithm is more favorable for large value of N, which is 
quite typical in OFDM systems. 

Algorithm 1 Proposed dual decomposition algorithm for (O 
1: Input the system parameters {TV, Pi ma x, p2max}, the chan- 
nel coefficients {h\ n , h 2n }^ =1 , and a solution accuracy e. 
Set fc = 1; Initialize dual variables u l = 1. 
repeat 

Compute the optimal {Pi„,P 2n } according to ( fT2b 
and 03) for V n G N\ 

Update the dual variables i/ k according to d2"Tk 
fc := fc + 1; 
until \\v k - <e||i/ fc - 1 ||. 

Output the optimal primal solution {Pf, P 2 } and PJ^a = 
min{C M A,i( P *, P 2), « = 1,2,3}. 



V. Simulation Results 

We consider an OFDM system with N = 32 subcarriers. 
The frequency-domain channels are generated using 8 inde- 
pendently and identically distributed Rayleigh distributed 
time-domain taps with unit variance |6|. The separate power 
constraints are set as Pi roax = -P 2 max = ^j?max, and ft = 0.5. 

Our proposed multi-subcarrier DF relay strategy is denoted 
as "Type 1 DF" scheme. Two reference schemes are considered 
in our simulations: The first one is the per-subcarrier two- 
way DF OFDM relay strategy in 0, which is denoted as 
"Type 2 DF" scheme; the second one is the two-way AF 
OFDM relaying scheme with optimized tone permutation in 
[6|. We divide the sum rate (approximated by the lower bound 
2i?x in Type 1/2 DF scheme) by N and use this per-subcarrier 
sum rate to evaluate performance at different average SNRs, 
which are only related with the power constraints Pi msK 's. 

Fig Elpresents the performance of different two-way OFDM 
relay strategies. The best performance is achieved by Type 1 
DF scheme with optimal power allocation (PA). At the spectral 
efficiency of 2 bits/s/Hz, Type 1 DF scheme with optimal PA 
provides a coding gain of about 2.5 dB compared with Type 2 
DF scheme, by performing channel coding across subcarriers. 
The PA gain between optimal PA and uniform PA of Type 1 
DF scheme is given by 1.6 dB. It is interesting that Type 1 DF 
scheme with uniform PA even outperforms Type 2 DF scheme 
with optimal PA, when the average SNR is in the region [0 dB, 
20 dB]. 

Although Type 1 DF scheme has no advantage over the 
AF scheme in the high SNR region due to its additional fully 
decoding requirement at Tr, it outperforms the AF scheme 
in the low and median SNR region. The intersection of the 
curves for Type 1 DF scheme and the AF scheme is at about 
17.5 dB, which is 5 dB higher than that for Type 2 DF scheme 
and the AF scheme. 

VI. Conclusion 

We have proposed a novel DF relay strategy for two-way 
OFDM relay networks and derived its achievable rate region. 
The key idea is making use of cross-subcarrier channel coding 
to fully exploit frequency selective fading. An efficient duality- 
based power allocation algorithm is also proposed to maximize 
the symmetric exchange data rate in both directions. Our 
simulation results suggest that the proposed DF strategy has 
better performance than existing DF or AF two-way OFDM 
relay strategies in the moderately low SNR region. We believe 
this two-way DF strategy tends to be optimal, i.e., achieving 
the capacity region outer bound, in the moderately low SNR 
region. The optimality of the proposed two-way DF strategy 
and the effect of channel uncertainty are currently under our 
investigation. 
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